Transport characteristics of nanojunctions far-from-equilibrium 
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We study the tunneling transport through a nanoj unction in the far-from-equilibrium regime 
at relatively low temperatures. We show that the current-voltage characteristics is significantly 
modified as compared to the usual quasi-equilibrium result by lifting the suppression due to the 
Coulomb blockade. These effects are important in realistic nanojunctions. We study the high- 
impedance case in detail to explain the underlying physics and construct a more realistic theoretical 
model for the case of a metallic junction taking into account dynamic Coulomb interaction. This 
dynamic screening further reduces the effect of the Coulomb blockade. 

PACS numbers: 73.63.-b, 73.63.Rt, 73.23.Hk 



I. INTRODUCTION 

Great efforts in contemporary materials science re- 
search focus on transport properties of advanced nano- 
materials. In particular in arrays of nano particles, the 
interest is motivated by the fact that these can be treated 
as artificial solids with programmable electronic proper- 
tieslil. The ease of adjusting electronic properties of gran- 
ular materials is one of their most attractive assets for 
fundamental studies of disordered solids and for targeted 
applications in nanotechnology. The parameters of gran- 
ular materials are in many ways determined by the prop- 
erties of their building blocks: grains and tunnel junc- 
tions. The equilibrium properties o£. single grains and 
single junctions are well understood!! However, much 
less is known about the far-from-equilibrium properties of 
those systems, by which we mean that the system proper- 
ties cannot be described by just a perturbed equilibrium 
(or quasi-equilibrium) considerations. The understand- 
ing of far-from-equilibrium effects in tunnel junctions, 
the building blocks of most advanced nano-materials, is 
especially important for practical applications. Indeed, 
recent experimental research has focused on instabilities 
in the current- voltage characteristics -showing clear devi- 
ations from quasi-equilibrium resultslm This defines an 
urgent quest for a quantitative description of far-from- 
equilibrium properties of a single tunnel junction. 

In this paper we investigate the far-from-equilibrium 
cur rent- volt age characteristics of a tunnel junction (see 
Fig. |l|a). Electron transport in tunnel junctions is en- 
sured by the energy exchange between the tunneling elec- 
trons and energy reservoirs: since the electronic energy 
levels at the leads are unequal, tunneling is only possible 
if a subsystem of excitations capable of accommodating 
this energy difference exists. At not very high temper- 
atures, where the phonon density is small, the role of 
the energy reservoir is played by an electromagnetic en- 
vironment comprised of electron-hole pairs self-generated 
by the tunneling electrons. Here we concentrate on this 
low-temperature situation where phonons (bath) are ir- 



relevant for the tunneling transport. In our approach 
the interaction time between electrons and environment 
needs to be much smaller than the one between envi- 
ronment and bath, in order to have a fully developed 
environment. 



II. PHYSICAL DESCRIPTION OF 
NON-EQUILIBRIUM EFFECTS IN A 
NANOJUNCTION 



We start by expressing the tunneling current through 
a single junction as the difference of the electrons go- 
ing from the left [L] to the right [R] electrode and the 
ones traversing the junction from right to leftBH (see also 
Fig. |a) 
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FIG. 1. (color online) a) Sketch of the tunnel junction with 
two leads and b) illustration of the electron-hole pair gener- 
ation across the tunnel junction. These virtual bosonic ex- 
citations together with the probability to interact with the 
environmental modes inside the junction determines the total 
tunneling current. 
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FIG. 2. (color online) Current- voltage characteristics of nano- 
junctions at low temperatures (T/E c = 0.1, with E c being the 
Coulomb energy of the nanoj unction) for the equilibrium case 
"eq." (red, dashed curve, see also |H|) and the non-equilibrium 
case "non-eq." , which shows a clear enhancement of the cur- 
rent and reduced suppression due to Coulomb blockade at 
low voltages. This is in particular clear in the plot of the 

ratio 5 = Inon-eq.,highim P .(K)//eq.(XO in tne inset - The lin- 
ear (green, solid) curve shows the high temperature ohmic 
regime, which is independent of non-eq. or eq. considerations. 
The two non-eq. curves correspond to the high impedance 
case (blue, solid) and the dynamic Coulomb interaction case 
(dashed, turquois); see text for detailed explanations. 
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where f £ L ' R ^ are the electronic distribution functions 
within the leads, Rt is the tunnel resistance, and V is 
the voltage difference across the junction. The terms 
and (1 — f^) correspond to the occupied electron state 
with energy e in the left lead and the hole state with en- 
ergy e' in the right lead, respectively. The function P(s) 
determines the probability that the tunneling electron 
loses [gains] the energy e to [from] environment modes in 
the junction. In general, both f^ L,R ^ and P(e) are out- 
of-equilibrium functions. In the second line of Eq. (||) 
we introduced a bosonic distribution function, N £ (eV), 
which describes an electron-hole excitations across the 
junction, Fig. ^p). Explicitly it is given by N £ (eV) = 

(s + eV)-' /r *» fi% +eV)/ 2 (l " e\s +e vy.) ■ K 
the distribution functions at the electrodes are Fermi 
functions with equal temperatures T, then N £ (eV) = 
A/#(£ + eV, T), with Nb(s, T) being the equilibrium Bose 
distribution function. 

The corresponding backward rate in Eq. (Q) is 



+ CO 

/ de(eV-e)N_ £ (eV)P(-e). Here we note, that the 

— CO 

bosonic distribution function N £ (eV) depends in general 
on both lead temperatures; and the bosonic form of the 
backward rate 1^( RL ) can be written in this form only if 
f £ L ^ and f £ R ^ have the same functional dependence on 
energy and temperature. 

Using the bosonic description of the environment, the 
current- voltage characteristics Eq. ([[]) can be written ex- 
plicitly as 

/ + CO 
de (e-eV) {N e (-eV)P(s) + N_ e {eV)P{-e)} . 
-CO 

(3) 

This expression is general and we only assumed 
that the electron distribution functions in the contact 
leads have the same functional dependence on tem- 
perature and energy. In the out-of-equilibrium sit- 
uation the tunneling electrons also interact with the 
environment inside the junction. This environment 
is self-generated by the tunneling electrons and ther- 
mal fluctuations, and its influence on the tunneling 
transport is implicitly taken into account through the 
probability function P(e) in Eq. (||). In the case of 
heat transport, the situation is different, since heat 
can dissipate inside the junction, see Ref. [Q and the 
environment needs to be taken into account explic- 
itly, having its own distribution function n env (s,eV) = 
[(e - eV)N e (-eV) + (e + eV)N £ (eV)\ /(2s), Refs. @|. 

As mentioned before, if the leads are in equilibrium 
(which is mostly the case due to their bulk nature) and 
in addition the temperatures at the leads are the same, 
N e (eV) in Eq. ([|), becomes a Bose distribution func- 
tion. In the following we assume that this is the case, 
but emphasize that the junction environment is always 
far-from-equilibrium if a finite voltage is applied. 

Besides the distribution function 7V e (eV), the main 
distinguishing feature of our out-of-equilibrium consid- 
eration is the presence of the probability function P(e), 
which takes into account the interaction with the envi- 
ronment, determining the probability for the tunneling 
electrons to exchange the energy £, which is the excess 
energy of the electrons compared to the potential dif- 
ference of the leads. It is clear that this probability 
should decay for large energies and have a maximum 
when the energy matches the energy at which the en- 
vironment resonates. A second observation is that at 
large voltages this probability should get smeared out 
and the current is mostly determined by the distribu- 
tion function N £ (eV) rather than P(e) which therefore 
determines only the resistance of the junction in the 
Ohmic regime. Furthermore, high temperatures also 
broaden the probability. Both effects can be conveniently 
described by the introduction of an effective electron 
temperature T e = (eV/2) coth(eV/(2T)) [derived from 
T e = lim e _»o ^env(£j eV)], which is equal to T for small 
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FIG. 4. (color online) Far- from- equilibrium probability func- 
tion P(s) at different voltages eV/E c — 0,0.5, 1 for both the 
high-impedance and the dynamic Coulomb interaction cases. 
The temperature is T/E c — 0.1 and a — 100, where the pa- 
rameter a = e 2 dv is defined below Eq. (^). Here E c is the 
Coulomb energy of the nanoj unction. 



FIG. 3. (color online) Far-from-equilibrium current- voltage 
characteristics of nanoj unction depending on voltage and 
junction temperature. Here E c is the Coulomb energy of 
the nanoj unction. One clearly sees the crossover from the 
suppressed current at small voltages and temperatures to the 
Ohmic regime at high temperatures. 



voltages and proportional to eV/2 for eV ^> TM There- 
fore, the probability function P(s) can be approximated 
by a Gaussian function where the electron temperature 
T e determines its width in the high resistive case, which 
we will discuss in detail in the next section. 



III. NANOJUNCTION WITH HIGH 
IMPEDANCE ENVIRONMENT 

We now turn to the experimentally important case of 
an environment with a high impedance as compared to 
the quantum resistance, Rq. In this limit, the tunneling 
electrons easily excite the environment modes. The prob- 
ability function P(e) for electron-hole pairs with energy 
e to appear in the junction in Eq. (|J) can be written as 



P{e) = (l/V^A 5 ) exp [-(s - 2E C ) 2 /2A 2 ] . (4) 

Here A = 2(E c T e ) 1 / 2 is the characteristic width of the 
distribution function with E c being the Coulomb energy 
of the nano tunnel junction. We note that this form of 
the probability function P(e) in Eq. (||) depends on the 
electron temperature T e and not on-jthe lead temperature 
T, as in the quasi-equilibrium caseH. 

Substituting this function P(e) into Eq. (^|), we ob- 
tain our first main result for the current-voltage char- 
acteristics of a tunnel junction. In particular, Fig. || 
represents the I — V characteristics at low temperatures 
(T/E c = 0.1) for the equilibrium case "eq." (red, dashed 



curve, see also Ref. g) and the non-equilibrium case "non- 
eq." . It shows a clear enhancement of the current and re- 
duced suppression due to Coulomb blockade at low volt- 
ages. This is in particular clear in the plot of the ra- 
tio S = / n on-eq.(^)/%q.(V) in the inset. The (green, 
solid) linear curve shows the high temperature ohmic 
regime, which is independent of non-eq. or eq. consider- 
ations. The full temperature and voltage dependence of 
the current-voltage characteristics of tunnel junction is 
shown in Fig. ||. This figure clearly shows the crossover 
from the suppressed current at small voltages (eV <C E c ) 
and temperatures (T <C E c ) to the ohmic regime at high 
temperatures. 



IV. NANOJUNCTION WITH DYNAMIC 
COULOMB INTERACTION 

Next, we discuss the current- volt age characteristics, 
Eq. (|lj), of a nano tunnel junction comprised of two thin 
two-dimensional (2D) disordered conductors (leads) tak- 
ing into account the effect of Coulomb interaction explic- 
itly. To this end we need to calculate the distribution 
function P(e), appearing in the tunneling rate, Eq. (||) 
from first principles. 

In general this function can be written as P(e) = 
dt exp[J(t) + iet], where the function exp[J(t)] ac- 
counts for the interaction with the Bosonic environment. 
The far-from-equilibrium function J(t) can be written 
aJT 

CO 

J M= [ ^ p(a; )[iV w e-* + (l + iV w ) e -^-B w ], 

Z J CU 

-1 

(5) 

where the terms proportional to N u and 1 + N u cor- 
respond to the absorbed and emitted environment ex- 




FIG. 5. (color online) Far- from- equilibrium phase function 
J(£), Eq. (|), at different voltages eV/E c = 0,0.5,1 for 
the dynamic Coulomb interaction case. The temperature is 
T/E c = 0.1 and a = 100, where the parameter a — e 2 dis is 
defined below Eq. (|§|) and E c is the Coulomb energy of the 
nanoj unction. The real and imaginary parts are plotted sep- 
arately and the inset shows the behavior of the real part of 
function J(t) for small times (in the zero voltage case). The 
latter demonstrates that the high-impedance expansion works 
in this regime. 



citations, respectively, and B u = 1 + 2N UJ . (Here we 
concentrate on the simplest case when the temperature 
of absorbed and emitted excitations is the effective elec- 
tron temperature T e determined by the environment.) 
In equilibrium reduces to the Bose-function and the 
functional P(oo) recovers the result of Ref. The en- 
ergy relaxation time r e in the expression for J(t) deter- 
mines the low energy cut-off, since the electrons start to 
equilibrate on larger time scales, i.e. the non-equilibrium 
description does not hold anymore. 

The spectral function p(uj) in Eq. (||) is the probability 
of the electron-environment interaction. We assume that 
leads are identical and have the same diffusion coefficients 
V {L) V {R) v and densities Q f states, z/ L ) 

v. For the dynamic Coulomb interaction the spectral 
function p(uo) can be found, following Ref. as 



(Vq 2 -iuj) 2 



where i,j = 1,2 are the lead indices for the left and 
right side respectively, and U^(q_,cj) are the dynamically 
screened Coulomb interactions within (across) the elec- 
trodes. The form of spectral probability p(uj) [p(u) = 
2pi2 + pn + P22] depends on the structure of the environ- 
mental excitations spectrum and, thus, on the external 
bias. 

The screened Coulomb interaction in Eq. ([|) in 
Fourier space has the form 0q,w) = {[£/ (0) (q,^)]" 1 + 
2(q, ^)} - \ where ?7^(q, uj) = u(q)£+v(q)a x is the bare 
Coulomb interaction and 2(q, uS) the polarization matrix 
respectively with V\\ = vV q 2 (V q 2 — ioj)~ 1 5[y 

Below we consider quasi-two-dimensional (2D) infinite 



FIG. 6. (color online) Spectral density function p(uS) for dif- 
ferent values of dimensionless parameter a = 1, 10, 50, 100 in 
double-log representation. The parameter a — e 2 dis is defined 
below Eq. (Js|) - The data are plotted for the case when Thou- 
less energy, Eth is equal to the Coulomb energy, Eth = E c . 



leads meaning that a < I <C L, where a is the lead thick- 
ness, I the electron mean free path, and L the lead size 
in the x and y directions. In this case the bare Coulomb 
interaction has the form 



with = (1/2 — 8i\)d and d being the nano junction 
size (distance of the contacts), leading to u(q) = 2ite 2 / q 
and v(q) = 2ite 2 e- qd /q. 

The dimensionless matrix elements Uij of the dynam- 
ically screened Coulomb interaction (in units of e 2 d) are 
then given by 
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dzj dz. 



I I 5 V V 



4n 



q x 2 (<l) - coth" 2 (g) ' I#J X(q) coth(q) 



(8) 



where q = dq and Co = uj(d 2 / V) with the dimensionless 
function x(o) = 1 + coth(g) + and a = e 2 dis. Using 

these notations, we can write Eq. @ as 



p(uj) 



2e 2 d 



U 



Im / qdq- 
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l-(x(q) coth(g))" 



(q 2 — iuo) 2 



(9) 

The spectral function p is plotted in Fig. g| as a func- 
tion of frequency for different values of the dimensionless 
parameter a = 1, 10, 50, 100. Notably, the p- function de- 
pends only weakly on frequencies in the low frequency 
limit and decays algebraically as p ~ 1/cj 1 / 2 at very high 
frequencies. Here, we remark that the parameter a has 
a typical value of 100 in the metallic case we are consid- 
ering here. 

Using Eq. (^) we numerically evaluate the behavior of 
the J(t)-function, Eq. (^|), which accounts for the inter- 
action with the environment. Its behavior for a typical 
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parameter value of a = 100 at different electron tem- 
peratures is presented in Fig. ||. The imaginary part of 
J-function is antisymmetric and almost voltage indepen- 
dent, while the slope of the real part is voltage depen- 
dent. The imaginary part contributes to the oscillatory 
factor in the expression for probability function P(e) in- 
troduced above Eq. (||). The behavior of real part is 
more important, since it describes the interaction with 
the environment and makes the P-integral convergent. 
At small dimensionless times (tE c <C 1) the real part of 
J(t) function has quadratic behavior (see the insert in 
Fig. ||), which corresponds to the high impedance limit 
for the environment discussed in detail before. However, 
at larger time scales the J-function shows linear behavior, 
which modifies the probability function P significantly. 

The behavior of P(e)- function is crucial for calculation 
of current- voltage characteristics in Eq. (|l|). We present 
the normalized probability function P(e) versus dimen- 
sionless energy in Fig. |j. In the case of dynamic Coulomb 
interaction, the P-function has two distinct features: i) 
a peak at low energies and ii) a long tail at high en- 
ergies. The first feature is related to the fact that in 
the limit of dynamic Coulomb interaction the screening 
effects are very pronounced and thus the original bare 
interaction E c is completely screened. This is in con- 
trast to the behavior of the P(£)-function in the high 
impedance environment, which has a peak at energies of 
order of Coulomb energy E c , see Fig. [I] for comparison. 
The second feature of probability function P(e), the ap- 
pearance of a long tail, increases the number of available 
states for energy absorption/emission of the environment 
enhancing the overall tunneling probability through the 
junction. This in combination with screening effect re- 
sults in a significant enhancement of the current at low 
voltages as compared to the high-impedance case for the 
metallic value of a = 100. If parameter a is decreased, 
i.e. the density of states lower, the Coulomb blockade 
gets restored, but also the resistance in the Ohmic regime 
increases. A current- volt age characteristics for a = 50 
is plotted in Fig. || (dashed, turquois), showing the en- 
hancement of the current at low voltages compared to 
the high impedance case. 

V. DISCUSSIONS 

Here we discuss the behavior of probability function 
P(e) in Eq. (||) and comment on the validity of our ap- 
proach at low temperatures. From Eq. (j|) follows that 
the probability function P(e) for zero temperature is 
proportional to the delta- function, P(e) ~ 5(e — 2E C ), 
meaning that no electron transport is possible below the 
Coulomb threshold. This is a consequences of our con- 
sideration of the tunneling transport in Eq. (H) being 



described in the lowest order in tunneling Hamiltonian. 
In this approximation higher order effects like electron 
co-tunnelingt2l is not taken into account. Co-tunneling, 
introduced in Ref. ||ll| , provides a conduction channel 
at low applied biases and temperatures, where otherwise 
the Coulomb blockade arising from electron-electron re- 
pulsion would suppress the current flow. The essence 
of a co-tunneling process is that an electron tunnels via 
virtual states thus bypassing the huge Coulomb barrier. 
There are two mechanisms of co-tunneling processes, 
elastic and inelastic. At very low temperatures only elas- 
tic co-tunneling exists meaning that electrons propagate 
through all virtual states without emitting/absorbing en- 
ergy. 

In this paper we only consider low (but not very low) 
temperatures were a bath (phonons) is inefficient (the 
typical validity temperature range would be between IK 
and 100K). Therefore our approach is valid when the 
interaction time between electrons and many-body exci- 
tations (environment) is much smaller than the one be- 
tween environment and bath, which is the case at not 
very high temperatures where the number of phonons 
is small. If the environment interacts strongly with 
the bath, relaxation is provided by phonons (bath) and 
P(e) = 5(e). In that case Eq. (|l|) reproduces Ohm's law. 

Last, we comment on two recent experiments of Refs. [|| 
and [| where instabilities in the current- voltage charac- 
teristics were observed. In these experiments the prop- 
erties of macroscopic systems, arrays of nano-j unctions, 
were studied and therefore our results can not be directly 
applied to these results. However, transport through a 
quantum nano- material can be reduced to the single junc- 
tion problem with an effective medium that plays also the 
role of a thermostat. This work is currently in progress, 
but requires in contrast to the singly unction also to con- 
sider heating effects of the mediumQ. 

In conclusion, we studied the tunneling transport 
through a nano-j unction in the far-from-equilibrium 
regime at relatively low temperatures. We showed that 
the current-voltage characteristics is significantly mod- 
ified as compared to the usual quasi-equilibrium result 
and demonstrated that for two cases: the high impedance 
and the dynamic Coulomb interaction case. One can ex- 
pect that our results will be important for electron trans- 
port in junction arrays, which will be the subject of a 
forthcoming work. 
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